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Abstract 

We report an exact analysis of a discrete form of the Chakrabarti- 
Stinchcombe model for earthquakes [Physica A 270, 27 (1999)] which 
considers a pair of dynamically overlapping finite generations of the Can- 
tor set as a prototype of geological faults. In this model the n-th genera- 
tion of the Cantor set shifts on its replica in discrete steps of the length of 
a line segment in that generation and periodic boundary conditions are as- 
sumed. We determine the general form of time sequences for the constant 
magnitude overlaps and hence obtain the complete time-series of overlaps 
by the superposition of these sequences for all overlap magnitudes. From 
the time-series we derive the exact frequency distribution of the overlap 
magnitudes. The corresponding probability distribution of the logarithm 
of overlap magnitudes for the n-th generation is found to assume the form 
of the binomial distribution for n Bernoulli trials with probability 1/3 for 
the success of each trial. For an arbitrary pair of consecutive overlaps in 
the time-series where the magnitude of the earlier overlap is known, we 
find that the magnitude of the later overlap can be determined with a 
definite probability; the conditional probability for each possible magni- 
tude of the later overlap follows the binomial distribution for k Bernoulli 
trials with probability 1/2 for the success of each trial and the number 
k is determined by the magnitude of the earlier overlap. Though this 
model does not produce the Gutenberg-Richter law for earthquakes, our 
results indicate that the fractal structure of faults admits a probabilistic 
prediction of earthquake magnitudes. 
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1 Introduction 



Earthquakes are outcomes of fault dynamics in the hthosphere. A geological 
fault is comprised of two rock surfaces in contact, created by a fracture in the 
rock layers. The two sides of the fault are in slow relative motion which causes 
the surfaces to slide. However, owing to friction the surfaces tend to stick and 
stress develops in the regions of contact. When the accumulated stress exceeds 
the resistance due to friction, the fault surfaces slip. The potential energy of the 
strain is thereby released, causing an earthquake. The slip is eventually stopped 
by friction and stress development resumes. Strain continues to develop till the 
fault surfaces slip again. This intermittent stick-slip process is the essential 
feature of fault dynamics. The overall distribution of earthquakes, including 
main shocks, foreshocks and aftershocks, is given by the Gutenberg- Richter 
law [1112]: 

logio Nr(>f > M)=a~b M (1) 

where, Nr(A^ > M) denotes the number (or, the frequency) of earthquakes 
of magnitudes A4 that are greater than a certain value M. The constant a 
represents the total number of earthquakes of all magnitudes: o = logj^Q Nr(A1 > 
0) and the value of the coefficient 6 w 1 is presumed to be universal. In an 
alternative form, the Gutenberg-Richter law is expressed as a relation for the 
number (or, the frequency) of earthquakes in which the energy released £ is 
greater than a certain value E: 

m{£ > E)^ E^^lf^ (2) 

where /3 « 3/2 is the coefficient in the energy-magnitude relation PIE]. 

One class of models for simulating earthquakes is based on the collective 
motion of an assembly of connected elements that are driven slowly, of which 
the block-spring model due to Burridge and Knopoff [H] is the prototype. The 
Burridge-Knopoff model and its variants [HI IZ] have the stick-slip dynamics 
necessary to produce earthquakes. The underlying principle in this class of 
models is self-organized criticality '5'. 

Another class of models for simulating earthquakes is based on overlapping 
fractals. These models are motivated by the observation that a fault surface is 
a fractal object ^ 1101 ITTl 1121 113j . Consequently a fault may be viewed as a pair 
of overlapping fractals. Fractional Brownian profiles have been commonly used 
as models of fault surfaces E] . In that case the dynamics of a fault is 

represented by one Brownian profile drifting on another and each intersection 
of the two profiles corresponds to an earthquake (16j . However the simplest 
possible model of a fault — from the fractal point of view — was proposed by 
Chakrabarti and Stinchcombe |17j . This model is a schematic representation of 
a fault by a pair of dynamically overlapping Cantor sets. It is not realistic but, 
as a system of overlapping fractals, it has the essential feature. Since the Cantor 
set is a fractal with a simple construction procedure, it allows us to study in 
detail the statistics of the overlap of one fractal object on another. 

In this paper we study the outcome of discrete dynamics in the Chakrabarti- 
Stinchcombe model. Our first aim is to construct the time-series of overlaps of 
the n-th generation of the regular Cantor set on its replica with periodic bound- 
ary conditions, where one set shifts on the other in discrete steps. We use a 
finite generation of the Cantor set since self-similarity in natural objects extends 
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only over a finite range. In this model the overlap magnitudes and their loga- 
rithms correspond to the energies released and the magnitudes of earthquakes 
respectively. Our second aim is to obtain the statistics of the overlap magni- 
tudes. From the time-series we derive the exact frequency distribution of the 
overlap magnitudes. For the n-th generation of the Cantor set the probability 
distribution of the logarithm of overlap magnitudes is found to assume the form 
of a binomial distribution for n Bernoulli trials and not an exponential law or 
a power law of the forms expressed in Eqs. Q and (0). As the Cantor set is 
approached through successive generations, the most probable overlap magni- 
tude is found to approach the cube-root of the maximum overlap magnitude. 
Our third aim is to determine whether it is possible to predict the occurrence 
of an overlap of a particular magnitude after a certain interval of time from an 
overlap that is known to have occurred. For this purpose we apply the theory 
of conditional probability. We consider the simplest case, that of two consec- 
utive overlaps in the time-series, where the magnitude of the earlier overlap is 
preassigned and using this we calculate the probabilities of all possible mag- 
nitudes of the later overlap. Each conditional probability is found to follow a 
binomial distribution for a certain number of Bernoulli trials that is determined 
by the magnitude of the earlier overlap. The asymptotic form of the most prob- 
able magnitude of the later overlap is found to be inversely proportional to the 
square-root of the magnitude of the earlier overlap. We use the expression for 
the conditional probability to study three important cases: where the magnitude 
of the later overlap is less than that of the earlier one, where it is the converse 
of the previous case and where the magnitude of the two overlaps are equal. On 
the basis of the fractal overlap analogy between a pair of fault surfaces and a 
pair of overlapping Cantor sets, our results indicate that, from the knowledge 
of earthquakes recorded in the past it is possible to predict, with a probability, 
the magnitude of earthquakes at least for a short term in the future. 



2 The model: Two dynamically overlapping Can- 
tor sets 

We consider a discrete form of the model of fault dynamics proposed by Chakrabarti 
and Stinchcombe |17| . The construction of the model requires two identical 
generations of the regular Cantor set. The procedure of constructing a regular 
Cantor set begins with a line segment of unit length, called the 'initiator'. This 
line segment is divided into three equal parts and the middle part is removed to 
obtain the first generation; that serves as the 'generator' of the Cantor set |18j . 
The procedure is repeated for each of the two line segments of the first genera- 
tion to obtain the second generation and so on. Therefore the n-th generation 
of the Cantor set is a finite set of 2" line segments, each of length 1/3". If this 
procedure is repeated an infinite number of times the remainder set of discrete 
points is called the regular Cantor set. In the rest of this paper the term Cantor 
set will always mean the regular Cantor set that we have just described. It is an 
exact self-similar fractal of dimension log3 2. The construction of the first few 
generations is shown in Fig. ^ Since it is practically impossible to generate the 
regular Cantor set exactly, we use the set of line segments obtained after a finite 
number of generations. It is reasonable to work with a finite generation of the 
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Cantor set since in all naturally occurring fractals the notion of self-similarity 
only applies between a lower and an upper cutoff length p[B|. 

The model is defined in terms of two dynamically overlapping fractals which 
are represented by the n-th generation of the Cantor set and its replica. The 
sequences of line segments in the two finite sets of the n-th generation are 
considered as a schematic representation of the fractal profiles of the two surfaces 
of the fault. The magnitude of the overlap of one fault surface on the other is 
measured in this model as the number of line segments of one set that overlap 
exactly on line segments of the other set and we denote it by Yn for the n- 
th generation sets. The motion of the fault surfaces is simulated by shifting 
one of the sets relative to the other in the following manner. We assume that 
initially the n-th generation of the Cantor set and its replica overlap completely, 
i.e., every line segment of one set overlaps exactly on the corresponding line 
segment of the other set. Then we shift one of the two finite sets relative to the 
other in one particular direction in uniform discrete steps of length 1/3", that 
is the length of a line segment in the n-th generation. The length of a step is 
chosen to be 1/3" in order to ensure that each line segment of one set either 
overlaps exactly on one of the other set or does not overlap at all. The overlap 
magnitude Yn after every step is given by the number of overlapping segments. 
Consequently it is convenient to define time t for this dynamical process as a 
discrete variable and measure its value as the number of uni-directional steps 
by which one set has shifted from its initial position relative to the other. The 
time-series of overlaps is given by the sequence {F„(t)}. Due to the discrete 
nature of the shifting process the magnitude of the overlap has discrete values. 
Further the structure of the n-th generation allows the number of overlapping 
segments to be in powers of 2 only: Yn = 2""'^, k = 0, . . . ,n. We shall write 
the magnitude of overlap as 2""*^ in order to indicate the generation n to which 
it belongs. We use periodic boundary conditions which produces a periodicity 
in the time-series. In Fig. |21we plot the time-series of overlap for the first few 
generations of the Cantor set. These plots show that the time-series has the 
appearance of a self-affine profile. 



3 Construction of the time-series 

In the following we construct the time-series of overlaps Yn {t) of the n-th gener- 
ation of the Cantor set on its replica by determining all of its structural features. 
We assume that periodic boundary conditions are assigned to the finite sets. The 
structural features of the time-series consist of a periodicity, a symmetry and a 
set of distinct temporal sequences of overlap magnitudes. The most important 
are the sequences of constant magnitude overlaps. The complete time-series is 
derived by superposing the constant overlap sequences of all possible overlap 
magnitudes. 

(i) Since the initial condition is chosen as the complete overlap of the two 
sets, the initial overlap is maximum and its magnitude is given by: 

y„(0) =maxy„ = 2". (3) 

(ii) Since the length of each line segment in the n-th generation of the Cantor 
set is 1/3" and the boundary conditions are periodic, the time-series of the 
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overlap for the generation n repeats itself after every 3" time-steps: 

Yjt) = Ynit + 3"). (4) 
Subsequently whenever we shall refer to the variable t and its values it will mean 

i = tmod3", (5) 

such that we shall study the time-series within the period < i < 3". 

(iii) Since the structure of all generations of the Cantor set is symmetric 
about its center and the boundary conditions are periodic, the time-series of 
the overlap in the period < i < 3" (i.e., in all periods of the kind /c3" < t < 
{k + 1)3", fc = 0, 1, 2, . . .) has a symmetric form: 

y„(f) = r„(3"-<), o<t<3". (6) 

Therefore, for every sequence of overlaps occurring at time-steps {<} within a pe- 
riod of the time-series, we have a complementary sequence of overlaps occurring 
at the time-steps {3" — t} within the same period. 

(iv) After the first 3"^^ time-steps of shifting one n-th generation set relative 
to the other, the overlapping region is similar to the generation n — 1 of the 
Cantor set. Consequently the next 3"~^ time-steps forms a period of the time- 
series of overlap of two sets of generation n — 1. Within this period of 3"^^ 
time-steps, after the first 3"~^ time-steps the overlapping region is similar to the 
generation n — 2 of the Cantor set and therefore the next 3"^^ time-steps form a 
period of the time-series of overlap of two sets of generation n — 2. This process 
occurs recursively and therefore a period of the time-series of the generation-rt 
overlap is a nested structure of the periods of the time-series of the overlaps 
of all the preceding generations, for example, after 3"~^ + 3"^^ + • • • + 3"^*^ 
time-steps the overlapping region is similar to generation rt — fc of the Cantor 
set, for which, following Eq.©, the magnitude of overlap is given by: 

Yn(j2T--^ =2^-^; fc = l,...,n. (7) 

The sequence of overlaps generated by Eq. (O and the complementary sequence 
obtained by the symmetry property of Eq. ((HJ, along with the initial condition 
(Eq. forms the skeleton of the entire time-series (illustrated in part (a) of 
Fig- El- The self-affine profile of the time-series observed in Fig. 13 is due to this 
property of the overlap. However it does not provide any more detail of the 
time-series. 

(v) If the generation-n of a Cantor set is shifted by a single time-step relative 
to another, the resulting overlap is just a unit line segment - owing to the 
periodic boundary conditions the last line segment of the former set overlaps on 
the first line segment of the latter set. This overlapping segment is similar to 
the initiator of the Cantor set (i.e., generation zero) with respect to the ratio 
1/3". Similarly, after a shift by three time-steps, the last two line segments 
of the former set overlaps on first two line segments of the latter set and the 
overlapping region is similar to the first generation of the Cantor set with respect 
to the ratio 1/3"^^. In general, owing to the periodic boundary conditions, 
a shift of 3*^ time-steps (0 < fc < n) from the initial position results in the 
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overlapping of the last 2*^ line segments of one set on the first 2'"' segments of 
the other set and the overlapping region is similar to the generation-A: of the 
Cantor set with respect to the ratio 1/3"^'^'. Therefore, beginning with a unit 
overlap resulting after a unit time-step, a sequence of shifts (time-steps) that are 
in geometric progression with common ratio = 3 produce a sequence of overlaps 
with magnitudes that are in geometric progression with common ratio = 2: 

y„(3^-)=2^ fc = 0,...,n. (8) 

This structural feature is illustrated in part (b) of Fig. |21 The same geometric 
progression of overlap magnitudes also occur for a geometric progression of 
time-steps beginning with a unit overlap after the fifth time-step: 

y„ (5 X 3'=) = 2*=; fc = 0, ...,n. (9) 

For each of the sequences given by Eqs. (jHl) and the complementary sequence 
is obtained by the symmetry property of Eq. . 

(vi) Since the n-th generation of the Cantor set contains 2" line segments 
that are arranged by the generator in a self-similar manner, the magnitude 
of the overlaps, beginning from the maximum, form a geometric progression 
of descending powers of 2. As the magnitude of the maximum overlap is 2" 
(observed at t — 0), the magnitude of the next largest overlap is 2"^^, i.e., 
a half of the maximum. A pair of nearest line segments form a doublet and 
the generation-n of the Cantor set has 2"~^ such doublets. Within a doublet, 
each of the two line segments are two steps away from the other. Therefore, 
if one of the sets is shifted from its initial position by two time-steps relative 
to the other, only one of the segments of every doublet of the former set will 
overlap on the other segment of the corresponding doublet of the latter set, thus 
resulting in an overlap of magnitude 2"^^. An overlap of 2"^^ also occurs if 
we consider quartets formed of pairs of nearest doublets and shift one set from 
its initial position by 2 x 3 time-steps relative to the other. Similarly, in the 
case of octets formed of pairs of nearest quartets, a shift of 2 x 3^ time-steps 
is required to produce an overlap of 2""""^. Considering pairs of blocks of 2'"^ 
nearest segments (ri < n — 1), an overlap of magnitude 2"~^ occurs for a shift 
of 2 X 3"'^ time-steps: 

y„(2 X 3''i) = 2"-^ n =0,...,n-l. (10) 

The complementary sequence is obtained as usual by the symmetry property of 
Eq. ©. An illustration of this constant overlap sequence is given in part (c) of 

Fig.ia 

Because of the self-similar structure of the Cantor set the line of argument 
used to derive Eq. (|10|l will work recursively for overlaps of all consecutive 
magnitudes. For example, the next overlap magnitude is 2"^^, i.e., a quarter 
of the maximum. For each time-step ti at which an overlap of 2""^ segments 
occur, there are two subsequences of overlaps of 2"~^ segments that are mutually 
symmetric with respect to ti; one of the subsequences precedes ti, the other 
succeeds ti. Therefore the sequence of t values at which an overlap of 2"~^ 
segments occurs is determined by the sum of two terms, one from each of two 
geometric progressions, one nested within the other: 

y„(2[3'^i ±3''^]) = 2"-2; (11) 
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ri = l,...,7i-l; 
r2 = 0, ...,ri-l. 

The first term belongs to the geometric progression that determines the sequence 
of t values appearing in equation Eq. (|10|) for the overlaps of 2"~^ segments, 
while the second term belongs to a geometric progression nested within the first. 
The symmetry property of Eq. (jSJ provides the complementary sequence. 

In general the sequence of time-step values at which an overlap of 2"^*^ 
segments {1 < k < n) occurs is determined by the sum of k terms, one from 
each of k geometric progressions, nested in succession: 

y„ (2 [S"^! ± 3''^ ± • • ■ ± 3"'"-' ± 3"'"]) = 2"-'=; (12) 

k = 1, . . . ,n ; 

ri = fc — 1, . . . , n — 1; 
r2 = fc - 2, . . . ,ri - 1; 

rk-i = 1, . . . ,rfe_2 - 1; 
rk = 0, . . . ,rfc_i - 1. 

For each value of k in the above equation there is a complementary sequence due 
to the symmetry property of Eq. ©. Eq. ((T^ along with Eqs. Q, l@J and © 
determine the entire time-series of overlaps of the n-th generation of the Cantor 
set on its replica. Assuming that the initial overlap is given by Eq. ©, the 
time-series is the superposition of the sequences of constant magnitude overlaps 
given by Eq. H12() for all possible overlap magnitudes. 

Now we consider the special case of unit overlaps, i.e., when the overlap is 
only on a unit segment. We shall write the magnitude of a unit overlap for the 
n-th generation of the Cantor set as 2"""' in order to indicate the generation 
index explicitly. A unit overlap occurs when k = n in. Eq. H12(l . The sequence 
of t values at which these occur is given by: 

y„ (2 [3"-i ± 3"-2 ± • • • ± 3^ ± 3°] ) = 2"-" = 1. (13) 

The above equation shows 2"^^ occurrences of the unit overlap. The same 
number of unit overlaps also occur in the complementary sequence obtained by 
using Eq. (jBJ. Therefore, in a period of the time-series for the n-th generation, 
there are altogether 2" unit overlaps. 

4 Analysis of the time-series 

The analysis of the complete time-series determined by Eqs. Q, Q, © and 
Eq. 112|) is carried out in two parts. First we systematically derive the exact 
number of overlaps of any magnitude y„ in a period of the time-series. Next we 
apply probability theory. The probability of occurrence of an overlap of mag- 
nitude Yn after any arbitrary time-step is obtained directly from its frequency. 
We then derive the conditional probability of the occurrence of an overlap of 
magnitude if it is given that the preceding overlap in the time-series has 
magnitude Yn- 
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4.1 Frequency of overlap magnitudes 

In the following Nr(yrt) denotes the frequency of overlaps of magnitude F„, i.e., 
the number of times an overlap of Yn segments occurs in a period of the time- 
series for n-tli generation of the Cantor set. From the sequences defined by 
Eqs. ©, l(Tn|) and l(TT|l and their complementary parts, we get respectively: 

Nr(2") = l, (14) 

Nr (2"-i) = 2n (15) 

and 

n-l 

Nr(2""2) 2 ^ 2ri = 2n(n- 1). (16) 

ri = l 

Similarly we can calculate the frequency of overlaps of all magnitudes. For 
example, 

n— 1 ri — 1 

Nr (2"-3) = 2 ^ 2 j]] 2 r2 

ri=2 r2 = l 

= ^n(n-l)(n-2), (17) 

n— 1 ri — 1 r2 — 1 

Nr(2"-) = 2:^25:25:2^3 

ri=3 r2=2 r3 = l 

= ^„(7i-l)(n-2)(n-3), (18) 



n — 1 ri — 1 r2 — 1 ^"3 — 1 

Nr(2"-) = 25:25:25:25:2^4 

n— 4 r2— 3 ^3— 2 r4 — 1 

4 

= — n(n- l)(n-2)(n-3)(n-4). (19) 
15 

In general, from Eq. H12|l the frequency of overlaps of magnitude Yn — 2"^*^ is 
given by: 

n-l ri-1 '•fc_2-l 

Nr(2"-'=) - 2 5: 2 5: ••• 2 5: 2rk-i 
= Cfc r7.(n - l)(n - 2) • • • (n - fc + 1) 

= ^^(^T^- ^^«) 

The value of the constant Ck is determined from the case of unit overlaps in the 
following way. In the above equation we keep the index k constant and choose 
the generation index n — k. As a result we get the frequency of unit overlaps 
for the k-th generation: 

Nr (2'=-'=) = Ck k\. (21) 
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(We have followed the notation explained before Eq. that the magnitude of 
the unit overlap for the fc-th generation of the Cantor set is written as 2''~'' in 
order to indicate the generation index.) On the other hand, we get the frequency 
of the unit overlap for the fc-th generation from the sequence defined in Eq. fl^ 
by replacing the index n by k: 

Nr (2'=-'=) = 2''. (22) 
Comparing Eqs. H21|l and (|22|l we get the value of the constant: 

= ^. (23) 

Now we have the complete formula for the frequency of overlaps of magnitude 
Yn — 2""'"' in the time-series for the n-th generation of the Cantor set: 

Nr(2"-'=) - 



2M " ), fc = 0,...,n (24) 



fc! {n-k)\ 

n 

where ^ ^ ^ denotes the binomial coefficient. The frequency distribution 

Nr(y„) obtained in Eq. (|24|l above for the overlap magnitudes F„ = 2"^'^, 
fc = 0, . . . , n, is exact. It is our central result and all subsequent results in this 
paper will be derived from it. It differs from the earlier claims ^|^, that the 
frequency distribution follows a power law; possible reason for this difference 
will be mentioned in the discussion. 

Using the binomial theorem, we have for the entire period: 

n ^ / \ 

^Nr(2"-'=)=5]2M n =3", (25) 

k=0 k=Q ^ ^ 

since there are 3" time-steps in a period of the time-series. From Eq. H24() we 
can calculate the cumulative frequency of overlaps with magnitudes that are 
greater than a preassigned value, say, 2""*^, by a partial sum of binomial terms: 

k k ^ \ 

Nr (y„ > 2"-'=) = ^ Nr (2"-'-) = E 2'' " ) ^ (26) 

which does not have a closed form. Fig. 0] shows the frequency distribution and 
cumulative distribution of overlap magnitudes for three consecutive generations. 
For large overlap magnitudes, which is realized for large values of n and k <^ n, 
Eq. H24(l reduces to the following asymptotic form: 

Nr(2"-^-) ~ (27) 
^ ' fc! 

This is evident from the expressions for fc — 0, . . . , 5 in Eqs. (|14|I - H19() . Using 
the notations Yn = 2'^~^ and maxF„ = 2", Eq. H27fl for the n-th generation 
may be written as: 

, \ l+log, n 

max y„ ^ ^ 



Yn ) (log2[(maxy„)/y„])! ^^^^ 
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which shows that the factor 1/Yn, which was the one obtained in ^7], is over- 
whelmed by the factor l/y,l°^^" in the expression for Nr(y„) for large n and 
small k. 



4.2 Prediction probability of overlap magnitudes 

We now analyze the time-series of overlap from the point of view of probability 
theory. The treatment by probability theory is necessary to determine whether 
the occurrence of an overlap of a certain magnitude can be predicted from 
the knowledge of the magnitude of an overlap that has already occurred. We 
continue to consider the case of the n-th generation of the Cantor set overlapping 
on its replica. Let Pr(y„) denote the probability that after any arbitrary time- 
step t we get an overlap of y„ segments. For the general case, Yn — 2""*^, 
A; = 0, . . . , n, it is given by: 



Prf2"-'=1 



Nr (2"-'^) 
ELoNr(2"-'=) 

3" V k 



n 



n-k / \ 3 



(29) 



The final expression in the above equation has the form of the binomial distri- 
bution for n Bernoulli trials [201 with probability 1/3 for the success of each trial 
and therefore Pr (F„) stands for the probability of the case where the number of 
successes is given by log2l^. Since Pr (2""'^) is maximum for k = [2(n + 1)/3J 
(and also for k = 2(n-fl)/3 — 1 when n-f 1 is a multiple of 3), the most probable 
overlap magnitude in general is given by: 

y„ = 2"-L2(n+l)/3j (30) 

where the floor function [xj is defined as the greatest integer less than or equal 

to X ini- 

Since maxy„ = 2", we get the following asymptotic relation for large 
values of the generation index n: 

f„ ^ (maxr„)'/^ (31) 

Next we consider the conditional probability that an overlap of magnitude 
occurs after the time-step t + 1 if it is known that an overlap of magnitude Yn 

has occurred after the previous time-step for any arbitrary t. The conditional 

probability is given by: 

where Nr (y„, Y^) is the number of ordered pairs of consecutive overlaps (y„, Y^) 
occurring in a period of the time-series for the n-th generation. It follows from 
Eq. (|12|l that overlaps of magnitude 2""'"' in the time-series are immediately 
succeeded (i.e., after the next time-step) by overlaps of magnitude 2'', < r < A: 
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only and never by an overlap of magnitude greater than 2*^ . For y„ = 2" and 
— r ~ 0, . . . ,k we have: 

Nr(2"-^2-)= ) ( J )■ (33) 

Therefore, from Eqs. (|23l, and we get: 

Pr(2n 2"-*^) 



J_//cy fc = 0, 
2^ r y ' r = 0, . . . , /c 



r \2 V 2 



(34) 



Here we find that the conditional probability follows the binomial distribution 
for k Bernoulli trials with probability 1/2 for the success of each trial. Eq. H34|l 
shows that the expression of Pr(2''|2"~'') is independent of the generation index 
n. Therefore, for fixed k and r, it has the same value for all generations n, 
provided that < r < A; < n. Since the conditional probability Pr(2''|2"^'^) is 
maximum when r — \ {k + 1)/2J (and also for r — (fc + l)/2 — 1 when fc + 1 is 
a multiple of 2), the most probable overlap magnitude to occur next to an 
overlap of magnitude Yn — 2""*^ is given by: 

f,; = 2L(*^-+i)/2J. (35) 

For large values of the generation index n and large k, k < n, from the above 
equation we get the following asymptotic relation: 



We consider three applications of Eq. H34|l assuming in each case that Yn = 
2""'^ and Y^ is the overlap magnitude occurring next to Yn in the time-series. 
First, we find that 



PiiK<Yn\Yn) = ^Pr(2n2"-'=) 



r=0 



n 



= 1 forO<A:<-. (37) 

This implies that, for k < n/2, an overlap of magnitude F„ = 2""*^ is always 
followed by an overlap of equal or less magnitude Yn- Consequently the case of 
an overlap of magnitude greater than that of the preceding one (i.e., F,' > Yn) 
appears only for k > n/2. Second, we have 



PriY:,>Yn\Yn) = J2 Pills'' I 2' 



r—n — k 

k 



r—7i — k 



^ E ( ' ) for^<fc<„. (38) 



Since the final expression in the above equation involves a partial sum of bino- 
mial coefficients, it does not have a closed form, except for the trivial case k = n 
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and a few other special cases. Therefore it must be calculated numerically for 
specific values of n and k. Third, if the magnitude of overlap after a certain 
time-step is Yn = 2^~^, the probability that an overlap of the same magnitude 
also occurs after the next time-step is given by: 



— k I rt^l — k\ 



Pr (y^-y„ I y„) = Pr(2"-'=|2' 



< fc < n/2 



1 



2M n-k 



k \ ./o^.^. • (39) 



n/2<k<n 



In this way we can derive the probability of various cases where Yj[ is specifically 
related to ¥„,. 



5 Discussion 



In this paper we have reported the construction and analysis of the time-series 
of overlaps for the discrete Chakrabarti-Stinchcombe model where the n-th gen- 
eration of the Cantor set shifts on its replica, with periodic boundary conditions, 
in discrete uniform steps. The n-th generation of the Cantor set is effectively a 
fractal for lengths between 1 and 1/3" (or, between 1 and 3", if the length of 
a line segment in the n-th generation is taken as the unit) and therefore it has 
no characteristic scale within this range; this is also reflected in the structure of 
the time-series within a period. However the frequencies and the corresponding 
probabilities of the binary logarithm of overlap magnitudes follow a binomial 
distribution, which has a characteristic scale given by the most probable overlap 
magnitude. The existence of a most probable value in the frequency distribu- 
tion creates the possibility of predicting overlap magnitudes in the time-series. 
We have shown the utility of this simple model in predicting an event from the 
knowledge of the preceding event. Though it is believed that earthquakes can- 
not be predicted |22| . the analysis presented in the previous section indicates 
that the overlapping fractal structure of faults ought to admit probabilistic pre- 
dictions of earthquake magnitudes. It is left to be explored how the probability 
of predicting an overlap magnitude in this model increases with the width of 
the preceding interval of time-steps within which all the overlap magnitudes are 
known. 

The Chakrabarti-Stinchcombe model, proposed for studying the mechanism 
of earthquakes, is by far the simplest model containing the rudiments of the over- 
lapping structure of a geological fault. The authors of the original model Tz] 
analysed, using renormalization group argument, the case of the Cantor set shift- 
ing continuously over its replica with open boundary conditions and reported 
that the density of the overlap magnitudes Y, in the infinite generation limit 
(n oo), was in the form of a continuous power law: p{Y) (x 1/Y. However 
the discrete form of the model defined for finite generations of the Cantor set 
with periodic boundary conditions - the subject of this paper ~ has a binomial 
distribution of the binary logarithm of overlap magnitudes (Eq. [^or Eg. I29|l . 
If, instead of periodic boundary conditions, open boundary conditions are used, 
the time-series will be of finite duration, consisting of 3" time-steps; with the 
initial condition given by the position of maximum overlap at i = 0, as stated 
in Eq. (j2Jl, there will be no overlapping line segments of the two n-th genera- 
tion sets for all odd values of t and the frequency of each overlap magnitude. 
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except the maximum, will be only half of the value given by Eq. H24|l . For large 
generation indices the asymptotic form of the binomial distribution is given by 
the continuous normal approximation |20| . but it never acquires the form of a 
power law. The difference in the results of Ref. |17j and this paper is apparently 
due to the difference in the nature of the shifting process which is continuous in 
the former and discrete in the latter (with the size of each discrete step being 
equal to the length of the line segment in the finite generation). The authors 
of Ref. [T7' also report that their result is valid for all types of Cantor sets and 
in fractals of higher dimensions, e.g., in the case of two overlapping Sierpinsky 
carpets. Though the discrete forms of these cases are yet to be studied, it has 
been reported in that for two spanning clusters at the percolation threshold 
on a square lattice - where each cluster is a random fractal embedded in two 
dimensions - the overlap magnitudes follow a normal distribution. 

Finally we consider a variant of the Chakrabarti-Stinchcombe model defined 
as the n-th generation of the Cantor set overlapping on its complement in the 
unit line segment; the complement of the n-th generation set is obtained by 
replacing each line segment (of length 1/3") in the latter by an empty segment 
and each empty segment by a line segment. The overlap magnitudes in this 
variant model are given by Yn — 2" — 2''', k — 0, . . . ,n. The time-series of 
overlaps can be directly derived from Eq. (|12|l by replacing each overlap mag- 
nitude 2""*^ in the time-series of the original model with an overlap magnitude 
2" — 2"^*^. Consequently the frequency distribution of overlap magnitudes is 

given by Nr (2" — 2*^) = 2"^'^ ^ ) ' ^^'^ variant appears to be more realis- 
tic than the original model since each of the two parts of a fractured rock is 
complementary to the other. However, in both cases the frequency of overlap 
magnitudes are found to follow binomial distributions. 
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Figure 1: The 'initiator' (n ~ 0) and the first four generations in the construction 
of the regular Cantor set. The first generation (n = 1) provides the 'genera- 
tor'. This figure illustrates the process of constructing successive generations 
described in the text. 
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Figure 2: Three important features of the time-series are illustrated for the first 
four generations (n = 1,2,3,4) of the Cantor set. The features within a period 
of each of the first three generations appear successively nested within a period 
of the fourth, (a) The skeleton of the time-series. It is the nested structure of 
the maximum overlaps for all the generations, (b) A sequence of overlaps whose 
magnitudes are in geometric progression with common ratio 2 which occur after 
the time-steps whose values arc in geometric progression with common ratio 3. 
(c) The sequence of overlaps of a constant magnitude that is half the maximum, 
shown for a period for all four generations. 
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Figure 3: A period of the time-series of overlap magnitudes for the n-th gener- 
ation of the Cantor set overlapping on its replica, for the first four generations 
{n = 1,2,3 and 4) according to Eqs. Q, ® and It shows that a 

period of the time-series for each generation is nested within a period of the 
time-series for the next generation. Owing to this recurrence, the sequence of 
overlap magnitudes in a period of the time-series for the Cantor set (71 = 00) 
forms a sclf-affinc profile. 
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Figure 4: (a) The frequency distribution of the overlap magnitudes Yn = 2""*^, 
k — 0, . . . ,n for three consecutive generations n — 11, 12, and 13 according to 
Eq. (|24() . The points have been joined by hues only to make out the asymmetric 
Gaussian appearance of the distributions, (b) The cumulative distribution of 
overlap magnitudes for the three consecutive generations shown in part (a), 
calculated by using Eq. (|26|1 . The cumulative distribution for each generation 
has been scaled by the period of the corresponding time-series. 
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